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Abstract: We obtain a non-relativistic diffeomorphism invariant string action as a 
special limit of the Nambu-Goto action in an FLRW background. We use this action 
to study non-relativistic string dynamics in an expanding universe and construct 
an analytic model describing the macroscopic properties of non-relativistic string 
networks. The non-relativistic constraint equations allow arbitrarily small string 
velocities and thus a 'frustrated' equation of state for non-interacting strings can 
be obtained without the need of a velocity damping mechanism. Assuming that 
colliding string segments reconnect by exchange of partners, non-relativistic string 
networks exhibit scaling behaviour, but with enhanced energy densities due to the 
smaller average string velocity. Non-relativistic string networks can be relevant in 
several contexts in condensed matter physics and cosmology. 
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1. Introduction 

Progress in solving string theory in various backgrounds can be done by considering 
sectors of the theory that decouple from the rest of the degrees of freedom in a 
suitable limit. Such decoupled sectors are characterized by having an altogether 
different asymptotic symmetry compared to that of the parent string theory. A 
well known example to such a truncation is the BMN |^ sector of string theory in 
AdSs X 5*^. Once a consistent sector is found, a complete worldsheet theory with the 
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appropriate symmetries can be written down without further reference of the parent 
theory. 

Non-relativistic string theory (see also ||^) in flat space is another consistent 
sector of string theory, whose world-sheet conformal field theory description has the 
appropriate Galilean symmetry . Non-relativistic superstrings and non-relativistic 
superbranes ||^, |^ are obtained as a certain decoupling limit of the full relativistic the- 
ory. The basic idea behind the decoupling limit is to take a particular non-relativistic 
limit in such a way that the light states satisfy a Galilean-invariant dispersion rela- 
tion, while the rest decouple. For the case of strings, this can be accomplished by 
considering wound strings in the presence of a background S-field and tuning the 
-B-field so that the energy coming from the i?-field cancels the tension of the string. 
In flat space, once kappa symmetry and diffeomorphism invariance are fixed, non- 
relativistic strings are described by a free field theory in flat space. In AdSs x 5*^ [0], 
the world-sheet theory reduces to a supersymmetric free field theory in AdS2. 

It is an interesting question whether similar non-relativistic string actions can 
be constructed in an expanding spacetime and if so, whether non-relativistic strings 
could play a cosmological role in the form of cosmic strings. The study of cosmic 
strings^ has been catalysed in the past few years mainly due to theoretical motiva- 
tions, in particular the realisation that they are generic in Supersymmetric Grand 
Unified Theory (SUSY GUT) models and brane inflation [|, |T^. The latter pos- 
sibility is of particular significance as it provides a potential observational window 
to superstring physics |Tl], [1^. Further, the fact that the Planck satellite and laser 



interferometers such as LISA and LIGO may be able to probe a significant part of 
cosmic string tensions relevant to these models \T^\, opens the possibility of detecting 
cosmic strings in the foreseeable future. 

One can think of situations in which ordinary cosmic strings could behave non- 
relativistically. Network simulations in a matter or radiation dominated universe []T4 | 



suggest that, at late times, string segments move relatively slowly and coherently on 
the largest scales, but also show evidence that small-scale-structure |]15|, |1^ which 
is largely responsible for damping energy from the network through the formation 
of minuscule loops, remains relativistic as Hubble damping is inefficient at scales 
much smaller than the horizon [|T^. However, the situation is different for strings 



in de Sitter spacetime, where Hubble damping can be very efficient rendering the 
strings essentially non-relativistic. This may be relevant for late time cosmology 
as observations |jl8|, suggest that the universe is already entering a de Sitter 
phase. Further, non-relativistic string networks have been considered as Solid Dark 
Matter (SDM) pO|, ^ and more recently as an alternative explanation of galactic 



^The dynamics of cosmic strings can be described by considering perturbations around a static 
solitonic string solution. Keeping all orders in the perturbations results in a relativistic effective 
string action, while keeping only up to quadratic order gives rise to the non-relativistic string action 
we will consider in this paper j^]. 
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rotation curves. 

It would thus be desirable to have an effective diffeomorphism invariant action^ 
describing the dynamics of non-relativistic strings in a cosmological context. On 
the other hand, the fact that one can construct a consistent worldsheet theory of 
non-relativistic strings at quantum level (in flat space) also motivates the study of 
fundamental non-relativistic strings in an expanding spacetime. In this paper, we 
point out that a non-relativistic diffeomorphism invariant action can be obtained in 
the case of a Friedmann-Lemaitre-Robertson- Walker (FLRW) spacetime as a limit 
of the relativistic Nambu-Goto action, and study the dynamics and cosmological 
evolution of non-relativistic strings. 

The structure of the paper is as follows. In section we review some basic re- 
sults about the Nambu-Goto action and the dynamics of relativistic strings in an 
expanding universe, which will be useful for comparison to the non-relativistic case. 
In section ^ we obtain a non-relativistic diffeomorphism invariant worldsheet action 
by taking a particular limit of the Nambu-Goto action in expanding spacetime. We 
move on in section ^ to study non-relativistic string dynamics as described by this 
action. The physical interpretation of non-relativistic strings as well as their possible 
coupling to cosmology-in particular the effective equation of state of an 'ideal gas' 
of non-interacting, non-relativistic strings-are discussed in section ^. The effect of 
string interactions is left for section ^ where macroscopic models for the cosmological 
evolution of both relativistic and non-relativistic string networks are discussed. In 
section |^ we solve numerically the non-relativistic network model for a wide range of 
parameters and compare the results to those of the relativistic model. In section ^ 
we discuss possible applications of non-relativistic strings in condensed matter and 
cosmological contexts. Finally, we have three appendices which describe an alter- 
native derivation of our non-relativistic string action as a semiclassical expansion 

around the vacuum solution (appendix |^, the Hamiltonian formulation of rela- 
tivistic and non-relativistic string dynamics (appendix P) and the construction of a 
spacetime energy- momentum tensor for the non-relativistic string (appendix jCj). 



2. Relativistic String in Expanding Spacetime 

Let us first consider a string moving in a D + 1 dimensional spacetime with metric 
Gmn {M, N = 0,1,2, D). Its world history is described by a two-dimensional 
spacetime surface, the string worldsheet = x'^la^), i = 0,1. The dynamics is 
governed by the Nambu-Goto action 

SNG = ~TnJv^d^a, (2.1) 

^Note that Ref. ||2^ considers an action applicable to a 'continuous medium' with internal struc- 
ture, which is invariant under limited reparametrisations preserving the worldlines of the constituent 
particles. Here we consider a diffeomorphism invariant non-relativistic action. 
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where Tr is the string tension and 7 is the determinant of the pullback of the back- 
ground metric on the worldsheet, = G MN{,x)dix'^^ djX^ . 

The equations of motion for the fields obtained from this action are given by 



0, 



(2.2) 



ArA7 ■ ^i-^ 

where TfJ^^ is the {D + l)-dimensional Christoffel symbol and V^x^^ the covariant 
Laplacian of the worldsheet fields x'^ . 

By varying the action with respect to the background metric Gmn we obtain a 
spacetime energy-momentum tensor 

1 



rpMN^yA-, 



^^'^diX^'djx'' (5(^+^)(/ - ■ (2.3) 



The rigid symmetries of (|2.1| ) are given by the Killing vectors of Gmn- The Nambu- 
Goto action is also invariant under 2D diffeomosrphisms of the worldsheet coordinates 
0"*. One can use this freedom to fix the gauge by imposing two conditions on x^'^{a^). 

Now consider string propagation in an expanding Universe described by a flat 
FLRW metric 

Gmn = a{x°Yr]MN (2.4) 

in conformal time x^ = ri. 

A convenient gauge choice in this case is the transverse temporal gauge given by: 



XX = 
r = 

The equations of motion (2^) become |24|: 

e = -2^ei^ 

a 

i + 2f (1 -i2)i = 



(2.5) 



where 



„/2 



'-7 



,/2 



1 



1/2 



(2.6) 



(2.7) 



The variable e is related to the canonical momentum associated to the field x°(r). 
Indeed, in the transverse gauge 7oi=7t-o- = 0, we have: 



Tr 

PO = 



-77 



drX^ 



-TRaixyex"" 



{21 



which, after imposing the temporal gauge condition x° = 1, becomes — TRa(x°)^e. 

In the transverse temporal gauge, the energy-momentum tensor ( ^^.3] ) of a rela- 
tivistic string in a FLRW background is |T^: 

T*^^(r/,y^) = ^^^TRjda{ex''x^ -e-^x'^^x'^) 5^^\y^ - x\a,r^)) , (2.9) 
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where, having integrated out 6{ri — x^{t)), I runs from 1 to D. 

To construct the string energy one projects T*^^ on a spatial hypersurface rj = 
const, with induced metric h and normal covectors um = ('2(?7),0), integrating over 
the D spatial coordinates: 

E(r/) = -j VhnMnNT^''d''y 

^/hT^d^y. (2.10) 



Thus, due to the foliation, the energy can be constructed from the 00 component of 
the energy momentum tensor. Since, \/h = a{r])^ , equation ( |2.10D becomes: 

E{r^)= jTla{r])''d^y = air])TnJeda. (2.11) 

Therefore, the energy of the relativistic string is simply the tension times the physical 
string length, taking into account relativistic length contraction. 



3. Non- Relativistic Limit of Nambu-Goto Action in FLRW 

Now consider a string, charged under a background antisymmetric 2-tensor field 
Bmn, propagating in D + 1 FLRW spacetime. The string couples to B through a 
topological Wess-Zumino term, so that the total action reads: 

S = Smg + Swz = -TRj,/^d^a + qjB*, (3.1) 

where q is the string charge and B* the puUback of B on the worldsheet. We consider 
a relativistic string aligned in the directions, its transverse coordinates being 

X" with a = 2,3,...D. 

The non-relativistic limit p|, ^ Q of this string consists of rescaling the longitu- 
dinal coordinates 

x^" ^ujx'', fx = 0,1 (3.2) 

and taking the limit u oo. This yields a divergent term, coming from Sng, which 
(in some geometries) can be canceled by an appropriate choice of a closed Bmn- If 
we assume that the string is wrapped on a spatial circle 



X 



^ ~ + 27iR (3.3) 



then the chosen Bqi cannot be set to zero by a gauge transformation. 

The above procedure generally works in flat spacetime (in fact one only needs 
the longitudinal part of the metric be flat 0), but for a curved background it is 
not guaranteed that there is a choice of closed B which cancels the diverging piece 
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of the action. Non-relativistic superstring actions have been obtained in the case of 
AdSs xS'^. 

We will now see that the non-relativistic limit can also be taken in the case of a 
FLRW background. We write the Lagrangian density of the Nambu-Goto piece as: 



Cng = -TRj-det[g,, + Gab{v)^^X-^,Xb] 



= -TR^-detg,,^det[6) + g^''Gab{v)dkX-d,Xb] , (3.4) 

where gij = G^^diX^'djX" and Gabiv) = aivY^ab- 

Then, assuming a power law expansion a(r]) = r/"/^ (for example a = 2 resp. 4 
in radiation resp. matter dominated era) we obtain the non-relativistic limit of Sng 
by the rescaling ( p.2|) , which implies 

a{7]) cj°/2a(r7) . (3.5) 

Expanding the Lagrangian density in powers of the parameter uj we then obtain: 

Cng = -Trlu" |^V-det(? + ^^/^d^g''GM^^X%X' + O | . (3.6) 

We can then rescale the string tension by 

Truo'^^T^ (3.7) 

and take the limit uj —>■ oo, yielding a finite and a divergent piece: 

Tn 



Crc, = -Yy~'^''^99''GMdiX^d,X'' {3.1 



-Cdiv = -Tquj^ ^/ -detg = -TQUj'^a{7]Y ^J -dei{r|^y^iX^'^jX'') . (3.9) 

The divergent piece can be canceled by choosing an appropriate closed B^y. Indeed 
if we choose^ B^^ = a{rj)'^e^y the Wess-Zumino part of the Lagrangian becomes: 

]^uj\quj'^)a{T^fe^yt'^d,x^diXr (3.10) 

This term precisely cancels the divergent piece ( p.9| ) if one tunes the rescaled charge 
(gcu") with the string tension Tq. We are thus left with the Non-Relativistic string 
action: 

SNR = -^j ^f^g''Gab{v)^^X%X'd^a . (3.11) 



This action can also be derived by a 'semiclassical approximation' |^ from the 
classical solution: 

r T , M = 

= < A(T , M = 1 (3.12) 
[ , M = aG (2,...,D) 

(see Appendix ^ for details). 

^The chosen Bfj,^ is closed. Working in zweibeins e'^ we have: dB — ^d[a{x^)'^£f^ye^ A e"] = 
d[a(i°)2e°Aei] = 2ade° Ae° Ae^ +a2(de° Ae^ +e° Ade^) = a'^{-w°^ Ae^ Ae'^) - a'^e° Aw^° Ae° = 
where we have used de + w A e = 0, Cartan's structure equation with zero torsion. 
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4. Non-Relativistic String Dynamics 



The action ( p.ll| ) is characterised by 2D diffeomorphism invariance with respect to 
the worldsheet coordinates a* and global Galilei invariance (modulo time translations 
due to time dependence of the metric) with respect to the transverse spacetime 
coordinates X". The canonical variables satisfy two primary constraints: 

P.e'^Xu^'" + I [^G^\x') + T,X"^X"Ga,{x')^ = (4.1) 
p^x'^ + P,X'" = 0, (4.2) 

where p^, Pa are the canonical momenta corresponding to and X". 

Varying the action with respect to the transverse and longitudinal fields X" and 
x^ respectively, one obtains the equations of motion: 

di {y/^g'^djX'') + ^/^g'^Vl^diX^djX'^ + ^^g^^diX^djx"^ (doGbc) G''^ = 

(4.3) 



-gg Om^^ On^^ — — 



(4.4) 



subject to the boundary condition (|3.3|) . For the metric (|2.4| ) the Christoffel symbols 
vanish and the transverse equations of motion (|4.3|) relate the covariant divergence 
of the transverse fields X"" to the time derivative of the transverse metric Gab- We 
can use the 2D reparametrisation invariance of the action to fix the gauge. For our 
discussion it will be convenient to work in the static gauge 

a;0 - r = , - A(T = , (4.5) 

identifying worldsheet and background times, while allowing for multiple windings 
of the non-relativistic string. Indeed, defining a G [0,27r), the periodicity condition 
( [4.5| ) requires that 

X = nR, (4.6) 

where n is the string winding number. 

After fixing the gauge, the physical degrees of freedom of the non-relativistic 
string are the transverse coordinates X" and the corresponding momenta Pa- The 
equation of motion ( f4.3|) becomes: 

X" = -2-X'^ + X-^X'"" , (4.7) 
a 

which is the wave equation with a cosmological damping term — 2^X'^. This equation 
(for A = 1) has been used by Vilenkin EH] to describe small perturbations around 
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a straight cosmic string, and was obtained by taking the hmit -C 1, X'"^ <^ 1 of 
the relativistic equations of motion in the static gauge. Here, there is also a winding 
number A. One might be tempted to say that, for a/a = 0, equation ( [4.7| ) imphes a 
wave propagation velocity 

vl = . (4.8) 

However, one should remember that the physical coordinates are not a, r but rather 
= Act, x° = r, so rewriting (|4.71 ) in terms of the physical variables we get (in the 
case d/a = 0) 

dloX^ = dl^X^ , (4.9) 

which describes a wave propagating at the velocity of light. The non-relativistic 
string allows the propagation of waves along the longitudinal directions with the 
speed of light. However, the transverse velocities are not restricted, in contrast to 
the case of the relativistic string. 
An 'energy' 

'^0 = ^1 da (^^G'^'ix') + ToX'^^X^Gabix')^ (4.10) 

can be obtained from the constraint (|4.1|) , which in the gauge ( [4.5|) becomes: 

V^ = ]^jda (aToX'^X^ + A-^ToX'^X'^) G,,(a;°) . (4.11) 

This can be interpreted as the sum of the kinetic and potential energies of transverse 
excitations along the string. The actual string energy, obtained by integrating the 
projection of the energy-momentum tensor on a constant x^ hypersurface is (see 
appendix 0): 

£;(x°) = a(x°)-^Po. (4.12) 
Since ^^-translation is not an isometry of (|2.4|) the Lagrangian is not time- 
translationally invariant and po is not conserved. In fact, its time evolution can 
be found from the longitudinal equations of motion ( [4.4| ). In the gauge ( ^4.5| ) the 
= component of ( [4.4[ ) becomes: 

i (X"X„ + A-=^X'"X:) • = A-2 (x'"X„)' + ^ (a-2x"^X^ - X"X,) . (4.13) 
Integrating we obtain: 

Po = aaXT^j da {\-^X'^X'^ - X"X^) 5ah , (4.14) 



where the boundary term gives no contribution due to the periodicity condition ( |3.3| ). 
Similarly, from the constraint ( [4.2| ) we define the momentum Vi along the string 

Vi = -\l PaX'-da = -To / XaX"'da (4.15) 



A. 

in the gauge ( |4.5| ). Translational invariance then dictates that Vi is conserved, as 
can be easily verified using the equations of motion. 
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5. Physical Interpretation and Cosmology 



5.1 The NR Particle vs NR String Limit 

The non-relativistic hmit is generaUy understood as a low velocity limit, which can 
be formally obtained by sending the speed of light c to infinity. This procedure works, 
at least in the case of the point particle, although there are some conceptual issues 
involved when taking limits of dimensionful constants like c [26, 27|. A safer route 
is to keep c constant and rescale the time coordinate by a dimensionless parameter, 
say uj, taking the limit uj oo. One can thus obtain a reparametrisation invariant, 
non-relativistic action for the point particle. The naive application of this to the 
case of the string fails^ but this problem was solved with the realisation (see also 
[0]) that in order to obtain a Galilei invariant string action one has to rescale both 
longitudinal coordinates, not the time coordinate only. In a sense, one can speak 
of a non-relativistic 'particle' limit, obtained by taking f <^ 1 and a non-relativistic 
limit for extended objects for which one has to rescale all worldvolume coordinates, 
as we did in section ^ for the case of the string. 

The rescaling of the longitudinal string direction corresponds to the assumption 
{dy/dxY <^ 1, which one makes when deriving the wave equation by applying New- 
ton's 2nd law on an infinitesimal string segment. In the rescaling prescription we 
followed, the waves move along the string at the speed of light as the string ten- 
sion equals the mass per unit length. One usually thinks of non-relativistic strings 
as 'violin-type' having a small tension compared to their mass per unit length and 
thus a subluminal 'sound speed' along the string. In this sense, the strings we con- 
sider here are 'hybrid', having a relativistic speed of propagation along the string, 
but transverse Galilei invariance. However, it is precisely this hybrid action (in flat 
space) which arises in the simplest Lorentz invariant field theories when one studies 
the low-energy dynamics of domain wall solutions. Strings with subluminal propaga- 
tion speeds (which would correspond to a differentiation between the string mass per 
unit length and the tension) can arise in more complicated models, which allow for 
spontaneously broken longitudinal Lorentz invariance through a current generation 



mechanism on the string worldsheet^ |]30|, |3l|. To obtain such string actions as a 



non-relativistic limit of the Nambu-Goto action, one would have to rescale each of 
the longitudinal coordinates by a different factor and take both factors to infinity 
while keeping their ratio constant. 

Note that, in order to ensure that the antisymmetric field B used to cancel the 
divergent piece of the action cannot be gauged away, the non-relativistic string had 
to wind a compact dimension, say a;^ ~ a;^ + 2ttR. In fact, the divergent piece of 
the action is a total derivative with respect to the worldsheet coordinates P, so if 



■^The string obtained in this limit has a fixed length and no physical oscillations (see Ref. [p8|). 
^See Ref. for a discussion of the relation between strings with broken longitudinal Lorentz 
invariance and Kaluza-Klein strings in one dimension higher. 
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the action ( p.ll ) was to be interpreted as an effective non-relativistic action, one 



could simply drop this term without requiring that the string is wound. However, if 
the non-relativistic string is to be interpreted as a fundamental object, consistency 
requires a non-trivial winding. In this case, there are two distinct scales, namely Tq 
of dimension mass-squared, which appears in the action and the mass scale 

m = 27mi?To, related to the geometry (through the compactification radius R) and 
the string winding number n. In fact, when quantising the non-relativistic string [^], 
one encounters again the necessity of winding, as the mass m is needed to define 
the energy states of the non-relativistic string spectrum. In the fiat case there are 
no physical states with zero winding number [Q. Also note that in deriving the 
non-relativistic string action (|3.11| ) we have defined the tension Tq by a rescaling of 



the relativistic string tension Tr, appearing in the Nambu-Goto action (see equation 
&)■■ 

To = cu"T^, (5.1) 

where the expansion exponent a is positive and u is taken to infinity. Interpreting 
To as the physically relevant quantity which is to be kept constant, the relativistic 
tension Tr goes to zero as u tends to infinity. 

Finally, we comment on the stability of the non-relativistic string. A closed non- 
relativistic string is more stable to breakage than its relativistic counterpart^. This 
is a consequence of the winding, which only allows a discrete number of potential 
'splitting points' along the string. From an astrophysical perspective, ordinary cos- 
mic string loops decay through gravitational radiation, which mainly couples to the 
kinetic energy of the fluctuations. In particular, the power in gravitational radiation 
scales with the sixth power of the root mean square (rms) velocity (see for example 



1^, |3^). Thus, if such non-relativistic strings were to play an astrophysical role, 
their decay rate would be power-law suppressed. For long stings, the main energy- 
loss mechanism is through string intercommutation, which removes string length from 
the long string network. This is also expected to be suppressed for non-relativistic 
strings as the interaction rates are proportional to the string velocities. We shall now 
consider the possibility of coupling non-relativistic strings to cosmology. 

5.2 Coupling to Gravity and Cosmology 

The non-relativistic action we have analysed describes the dynamics of the indepen- 
dent degrees of freedom of the non-relativistic string, namely the transverse excita- 
tions. In obtaining this action we have introduced a closed B field, which cancels 
the divergent piece corresponding to the rest energy of the string. Alternatively, if 
we are not interested in quantisation, we can simply drop the divergent part of the 
action-without introducing the B field-because it is a total derivative (cf the case of 
the point-particle). Here, we will follow the latter approach. The energy-momentum 



^We thank F. Passerini for discussions of this point. 
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tensor of the non-relativistic string (Appendix 0) therefore describes the energy of 
the transverse excitations but does not include a contribution from the rest mass of 
the string. However, when one couples non-relativistic matter to General Relativ- 
ity it is necessary to include the rest mass mgc^ in the energy-momentum tensor, 
which gives the main contribution to the part while kinetic contributions are 
sub dominant. Following this logic we will add the rest mass of the string to the 
part of the energy-momentum tensor of Appendix |C|, which can then be coupled to 
Einstein's equations. From now on we work in D = 3 spatial dimensions. 

Consider a cosmological setup where the cosmic fluid has a component due to a 
gas of non-interacting, non-relativistic strings. To obtain the energy density of the 
string fluid, one has to sum the contributions of all string segments in the network 
and, as we discussed, it is the rest energy of the segments which will give the dominant 
contribution. This is in analogy to a gas of non-relativistic particles (dust), where 
the dominant contribution to the energy density is p = Tq = mQ-n + O(f^), where 
mo is the particle rest mass, n the rest frame number density and v the rms particle 
velocity. The off-diagonal terms of the energy-momentum tensor of the particle 
fluid average out to zero by summing over all particles with random velocities in all 
directions, whereas the T\ = —p components are proportional to the kinetic energy 
density, which for non-relativistic particles is negligible so that p <C p. 

In the case of a 'string gas' one can obtain an effective energy-momentum tensor 
in an analogous manner, by approximating the string network as a collection of 
straight string segments moving with average velocity v, and averaging over string 
orientations and directions of motion. Let us first consider the relativistic case. The 
effective energy-momentum tensor can be constructed by considering a straight string 
oriented in the z direction say, and Lorentz-boosting its energy-momentum tensor 
in the ±x and ±y directions, then averaging and repeating the same procedure for 



strings oriented in the x and y directions |Q. The result is: 



3L2 



/ 

(l-t;V) 

(l-t;V) 

\ (1 - ^2^2) J 



(5.2) 



where /i is the string tension, L the average separation between nearby strings in the 
network and 7 = (1 — t;^)^^/^ a Lorentz factor corresponding to v. From ( ^.2|) the 
equation of state can be read: 

-P = (Tl) = ^(7-^ - v'){n) = - 2v'){n) ^p= ~{1 - 2v')p . (5.3) 

A similar procedure can be followed for non-relativistic strings, which are gener- 
ally expected to have much smaller string velocities. Indeed, for relativistic strings 
the constraint x"^ + x'^ = in the conformal gauge imposes that critical points on the 
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string move with the speed of hght, but for non-relativistic strings the physical string 
velocities can take any value. One can thus obtain the equation of state for such a 
non-relativistic string gas by using Lorentz boosts with 7 = 1 or, alternatively, by 
performing transverse Galilean boosts instead. The result is again p = — i(l — 2v'^)p, 
but with the difference that one can safely assume f <^ 1, unlike the relativistic 
network case, where the strings oscillate relativistically at small scales, while there 
is no known mechanism which is efficient enough to damp these excitations. Indeed, 
Hubble damping is inefficient at scales much smaller than the horizon, and for large 
scales, of order the string correlation length, numerical simulations (see for exam- 
ple [Q) demonstrate that string segments move more slowly and coherently, but 
at speeds large enough to produce significant deviations from the equation of state 
w = p/p = —1/3. 

Note that one can apply an analogous procedure for strings which have a tension 
T smaller than their mass per unit length p {T < p). In this case the resulting 
equation of state is: 

P = -^[T//i(l - v') - v']p = -^-[vl - (1 + viy]p, (5.4) 

where we have defined the 'sound speed' along the string vq = ^JT j p. Equation (|5.4| ) 
can in general lead to both positive or negative equation of state with p > — p/3. 
This is in contrast to vacuum (non- interacting) cosmic strings with p = T, where the 
rms speed does not exceed 1 / \/2 so the equation of state is nonpositive ( |5.3|) with 
p > —p/3. However, this is to be expected because in the limit T the 'string' 
describes a line-like structure of dust particles with < p <^ p. In fact, taking t>o ^ 
in equation ( |5.4| ) gives p = pv'^/3, or, in terms of the kinetic energy density pk, 

2 

P = gPfc , (5.5) 

which is precisely the equation of state for a gas of non-relativistic particles, following 
from ordinary kinetic theory considerations. In connection to the discussion of the 
previous sections, obtaining this kind of non-relativistic string from the Nambu-Goto 
action involves a rescaling of the longitudinal directions by different factors, the ratio 
of which determines the propagation speed vq- 

Finally, note that this discussion only applies to a 'perfect' gas of non-interacting 
strings. String intercommutations typically result in the removal of energy from the 
network in the form of closed string loops, significantly altering the above picture. 
Thus, a frustrated string network, with w ~ — 1/3, p cx eventually dominates 
over matter or radiation, but turning on string interactions will result to a differ- 
ent equation of state. For abelian string networks, where interactions are efficient, 
the resulting scaling law is p oc t^^, where t is cosmic time, which scales like ra- 
diation in the radiation era and like matter in the matter era. The cosmological 
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evolution of non-relativistic string networks, including the possible effects of string 
intercommutation will be discussed in the next section. 



6. Velocity Dependent One-Scale (VOS) Models 

In this section we discuss analytic models for the evolution of macroscopic variables 
describing the large-scale properties of a string network. We will first review results 
for relativistic strings and then construct a macroscopic evolution model for non- 
relativistic strings, based on the action To set up the physical picture we 
briefly summarise Kibble's one-scale model which captures the basic qualitative 
features of network evolution. 

Monte-Carlo simulations of cosmic string formation suggest that to a good ap- 
proximation the strings have the shapes of random walks at the time of formation 
[p5[] . Such 'Brownian' strings can be described by a characteristic length L, which 
determines both the typical radius of curvature of strings and the typical distance 
between nearby string segments in the network. On average there is a string segment 
of length L in each volume and thus the density of the cosmic string network at 
formation is 

pL fi 

P = JJ = J2^ (6.1) 
where /i is the string mass per unit length, which for relativistic strings is equal 
the 'tension' Tr appearing in the Lagrangian. Assuming that the strings are simply 
stretched by the cosmological expansion we have p oc a{t)~^. This decays slower 
than both matter (oc a~^) and radiation (oc a~^) energy densities and so such non- 
interacting strings would soon dominate the universe. 

Now consider the effect of string interactions. As the network evolves, the strings 
collide or curl back on themselves creating small loops, which oscillate and radiatively 
decay. Via these interactions enough energy is lost from the network to ensure that 
string domination does not actually take place. Each string segment travels on 
average a distance L before encountering another nearby segment in a volume L^. 
Assuming relativistic motion [v ^ 1) and that the produced loops have an average 
size L, the corresponding energy loss is given by pioops ~ L~^pL. The energy loss 
rate equation is therefore 

dp , , 

6.2 

a L 

Equation ( |6.2D has an attractor 'scaling' solution in which the characteristic length 



L stays constant relative to the horizon dn ^ t |3J]. The approach of string networks 



to a scaling regime has been verified by high- resolution simulations |36, 14 



Equation (pT^ ) was derived on physical grounds and it only captures the basic 
processes involved in string evolution, namely the stretching and intercommuting of 
strings. It does not take into account other effects like the redshifting of string veloc- 
ities due to Hubble expansion. In fact, it neglects completely the evolution of string 
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velocities, making the crude approximation that they remain constant throughout 
cosmic history. However, we can construct a more accurate Velocity-dependent One- 
Scale (VOS) model, based on the Nambu-Goto action (|2.1|). 



6.1 Relativistic Strings 

The relativistic VOS model P7|, extends Kibble's one-scale model, abandoning 



the constant string velocity approximation and introducing an extra variable, the 
rms velocity of string segments, whose dynamics is governed-as we will see-by a 
macroscopic version of the relativistic equations of motion (|2.6| ). Although the sim- 
ple one-scale model captures most of the qualitative features of macroscopic string 
evolution, this correction is crucial for quantitative modelling. Indeed, the average 
string velocity enters linearly in the loop production term, which provides the main 
energy loss mechanism of the string network, and so the evolution of string velocities 
significantly affects the string energy density. The resulting VOS model is still very 
simple depending on only one free parameter^ but, remarkably, it has been shown to 
accurately fit numerical simulation data throughout cosmic history We briefly 



sketch how the model is constructed from the microscopic equations of section ^ 
This will be useful for comparison to the non-relativistic case. 

Consider the relativistic string energy defined in section |] (equation (|2.11| )): 



E{r]) = a{r])TR jeda 



and take the first derivative with respect to conformal time 77. Using the equation of 
motion ( p.6|) for e, one finds 

E = -(l-2v^)E, (6.3) 
a 

where f ^ = / ex^ da/ Jedcr= (x^) is the worldsheet average of the square of transverse 
velocities. For a network of strings the energy density p is related to the total string 
energy hj E cc pa{ri)^. Therefore: 

p E d d / r,^ , . 

- = --3- = -2- (l + t;2) . 6.4 
p E a ^ ^ ^ 



To this we add a phenomenological term p4|, p!7| describing the production of loops 



when strings collide and curl back on themselves. The resulting network density 
evolution equation is: 

p = -2-{l + v^)p-c^, (6.5) 
a L 



^Strictly speaking there are two parameters in the VOS model, the loop production efhciency 
c and the momentum parameter k. For the second parameter, however, there exists a physically 



motivated ansatz (6.15), which expresses it in terms of the rms velocity v{t). Once this choice is 



made, one is only left with the freedom of tuning c when trying to fit numerical simulations. 
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where c is the loop production efficiency related to the integral of an appropriate 



loop production function over all relevant loop sizes |Tj]. This is treated as a free 
parameter which can be determined by comparison to numerical simulations. 

In the VOS model, the rms velocity v appearing in equation (|6.5| ) is promoted to 
a dynamical variable whose evolution is given by a macroscopic version of the Nambu- 
Goto equation of motion (p.6|). This equation can be obtained by differentiating 
and eliminating x using the equation of motion. This introduces the second spatial 
derivative x" which corresponds to string curvature and can be expressed in terms 
of the mean curvature radius of the network. Differentiating f ^ and using equation 
we find: 

2vv = (x^y = 2(x ■ x) - 2- U-kY - (i^)) . (6.6) 

a 

The second term is of purely statistical nature and has the effect of 'renormalising' 
the coefficient of the -v^ term which will find later. It has been demonstrated 

a 

numerically [0 to have small magnitude and thus can be neglected. 

Keeping only the ffist term and using the equation of motion for x we find: 

rx-x"e^^(i(T fx ■ x'(e^^)' (icr a / , , a,\ , ^ 

Jeda ^—[^ 2-((xV{x*». (6.7) 

The second term vanishes due to the gauge condition x ■ x' = 0. Further, within our 
approximations (x^) ~ (x^)^ so the third term becomes 2^f^(l — v"^). For the ffist 
term we need to express x" in terms of the local curvature vector. We define 

ds = V^da = Vl - i^eda (6.8) 

and the physical (local) radius of curvature by 

— = (6 9) 

ds^ 7^ ' ^ ' 

where u is a unit vector. Then: 

da"^ ds"^ ds 

Due to the constraint x ■ x' = the second term vanishes on dotting with x so we 
have: 



±-x"e-'da= J±-^{l-±'')eda = a{ri){{±-u){l-±^)/n) Jeda. (6.11) 



We define the momentum parameter k ||3^ by the equation: 



((i■u)(l-i^)/7^) = ^(l-t;2), (6.12) 
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where TZ is now the average string radius of curvature, numerically close to the 



correlation length L for Brownian networks [ll7| , p^ , |39| . With this definition, equation 
( |6.7| ) becomes: 

y = ^kil - v^) - 2%il - v') . (6.13) 
iZ a 

Changing to cosmic time t, with dt = adr] and ' = a-^ we finally obtain: 

"^'^ 'l-v^)(^- 2Hv) , (6.14) 



dt ' ' \n 

where H = is the Hubble parameter. Note that, since 



dt 

\dr]J 




and the physical coordinates Xphys are given in terms of the comoving ones x by 
Xphys = ox, the rms velocity v has the interpretation of physical peculiar velocity of 
string segments. Equation (|6.14|) has therefore a clear physical meaning: the rms 



peculiar velocities of string segments are produced by string curvature and damped 
by cosmological expansion. 

The momentum parameter is a measure of the angle between the curvature 
vector and the velocity of string segments and thus it is related to the smoothness of 
the strings. As v increases towards relativistic values the accumulation of small-scale 
structure renders the strings wiggly. Velocities become uncorrelated to curvature 
and k decreases. In particular it can be shown analytically that for flat space, where 
t"^ = 1/2, the momentum parameter vanishes for a wide range of known solutions 

An accurate ansatz for the momentum parameter k for relativistic strings has 
been proposed in 



(6.15) 

satisfying ^(1/^2) = 0. 

Note that the fact that v = 1/V2 in flat spacetime, can be shown analytically 
for closed loops only, but for long strings it is observed in numerical simulations [|T^] . 
For expanding or contracting spacetimes, v is less or greater than 1 / a/2 respectively. 
Hence for an expanding universe, string velocities are subject to the constraint: 

t;2<i. (6.16) 

In a matter or radiation dominated universe, Hubble expansion is too weak to signif- 
icantly reduce string velocities, which remain close to 1/2 at short scales This 
limitation does not apply to non-relativistic strings. 
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6.2 Non-Relativistic Strings 

For the non-relativistic string tlie energy of the excitations is given by (see Appendix 
0): 

^exc = a{r^)]^ J da (/iX^ + /xA-^X'^) = a{r])-' Vo , (6.17) 

where X are the transverse string coordinates and we have defined the tension /i = 
XTq. To that we must add the string mass 

Eq = a{ri)fi J da , (6.18) 

so that the total energy is: 

E = Eo + Eexc = a{ri)fx j da + a{ri)-' Vo . (6.19) 
Then, differentiating with respect to conformal time ( ' = J^), we have: 

E = -Eo + {a-^Vo)- = -Eo - -E,^, + a'' Vq 
a a a 

^(i + -W'--vAeo, (6.20) 



a \ 2 2 ^ 

where we have used equations ([4.11|) , ([4.13|) and defined the rms quantities: 

V' = ^ (X^) , (6.21) 

and 

= ^^^^^ = (A-^X'2) = {{d..lLf) , (6.22) 

corresponding to the average velocity of string segments and the average magnitude 
of string tangent vectors. The latter quantity parametrises small-scale perturbations 
on the string, W = corresponding to strings which are straight at scales smaller 
than the correlation length®. Thus, the term in equation ( |6.20| ) corresponds 

to the average elastic energy due to short-scale string deformations. In the non- 
relativistic limit one has W'^ <^ 1. 

Defining the energy density p oc Ea~^, and using 

E E p a , , 



we find 



a 



p = -- (2- + ^V'' ] p-cV^ , (6.24) 



*With this interpretation, one expects that W should have the effect of reducing the effective 
radius of curvature of the network. As we will see later, this is indeed the case. 
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where we have included a phenomenological loop production term, as in the rela- 
tivistic case. 

From ( |6.21|) we have 

2VV = (X^)- = 2(X ■ X) - 2^ - (X')) (6.25) 

as before. We neglect the statistical terms and using the non-relativistic equation of 
motion ( [4.7| ) we find: 

VV = ^-^^ = ^- p- (6.26) 

j da j da a 

In order to express X" in terms of the string curvature vector we define: 

ds = + (9^iX)2da;^ = Wl + X-^X'^da (6.27) 

and the physical radius of curvature: 

d^Y _ai7]) 

where we have introduced the 3- vector Y = ('x^,X) and a unit 3- vector u. Now: 



X-- = ^ ^ A' (1 + A-^X") + AX- (6.2<,) 

da^ ds^ ds 

In this case, the second term will not cancel on dotting with X, because X ■ X' 7^ 
for the non-relativistic string. Instead we have two terms: 



Jx-X"da = Jx-^{1 + \-'X'^) da + \-' jx- X' (in Vl + A'^X'^) ' da . 

(6.30) 

For the first term we note that, since X is normal to (x^, 0) in Cartesian coordinates, 

. d^X ■ d'Y 

and so we can use equation ( |S.28| ) to write: 

yx . ^ (1 + A-2X'2) da = a{r])^il + W^) j da . (6.32) 

Here, in analogy to the relativistic case, we have defined a momentum parameter k 
by: 

'(1 + A-2X'2) (X ■ u)/7^^ = ^(1 + W^) . (6.33) 
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For the second term in ( |6.3C ) we have: 

■ X' (in Vl + A-2X'2) 'da = X-'- J±- X^ f^'^'l^^^ 
A-^ /fx . X') (X' ■u)^da + X-' [(X- X') X" da 
a{r]) J / da + 0{VW^) , (6.34) 



where we have used equation ( |6.29| ) and defined the parameter k' by: 

A-(x.X')(X'.u)/7^) = ^ (6.35) 

Putting all the terms together, equation (|6.26|) can be rewritten (in terms of cosmic 
time t) as: 

^ = ^{k + k"W')-2HV, (6.36) 

with 

k" = k + k'. (6.37) 

Equations (|6.24| ), (|6.36|) form the Non-Relativistic Velocity dependent One-Scale 
(NRVOS) model. In principle one should consider W as a third dynamical variable 
and try to derive an evolution equation, as in the case of V. As a first approxima- 
tion we will assume that time variations in W do not have a significant impact, W 
remaining always small, and we will treat it as a constant parameter. This approxi- 
mation will be tested in the next section, where we will solve the NRVOS equations 
numerically, for different choices of the W parameter. 

Finally, one comment is in order regarding the magnitude of the parameter k'. 
From its definition in equation ( |6.34|) one expects k' <^ k. Indeed, k' measures 
the average value of (X • X')(X' ■ u) the first factor of which contains uncorrelated 
vectors, while for the second factor string tangents will generally be normal to the 
local curvature vector. On the other hand k corresponds to the average value of X ■ u 
and these two vectors are correlated, at least for smooth strings/small excitation 
velocities. Given that the k" term in (|6.36| ) is already suppressed by a factor OiW"^) 
it is a good approximation to set k" c::^ k. Then, W has the effect of 'renormalising' 
the effective radius of curvature TZ '^/(l + W^) (or equivalently the momentum 
parameter k — ^ k{l + W^)), as may be expected from its interpretation as a short- 
scale structure parameter. 



7. Relativist ic vs Non-Relativistic Network Evolution 

In this section we solve numerically the NRVOS equations for a non-relativistic string 
network and compare to the relativistic case. The naive expectation is that non- 
relativistic networks are denser than their relativistic counterparts because the small 
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string velocities reduce the effect of the loop production term. Physically, the trans- 
verse excitations on strings are non-relativistic so fewer loops are produced per unit 
time due to string self-intersections. Long string segment interactions are also sup- 
pressed due to the low collision rate corresponding to small velocities. 

To close the NRVOS equations we need to specify an ansatz for the non-relativistic 
momentum parameter k. For a velocity dependent model like the one we have de- 
veloped, it is not consistent to treat A; as a constant parameter. Further, in the 
relativistic case, its dependence on the rms velocity v (equation ( |6.15| )) is important 



in determining the scaling values of the network variables. The functional depen- 
dence of the momentum parameter on v can be obtained by considering 'curvature' 
and 'bulk' contributions to string velocities, as explained in Ref. [^. Following the 
discussion in that reference we take: 

kiv) = koil-v'), (7.1) 

where fco is a constant. This has the same functional dependence as the low- velocity 



limit of k{v) in Ref. but here we have left the overall normalisation ko as a 
free parameter. This reflects the fact that the non-relativistic string limit is not 
merely a low-velocity one. There is a difference between slowly moving, straight, 
relativistic strings and wiggly, non-relativistic strings. The defining property of the 
non-relativistic string is that its transverse excitations be Galilei, as opposed to 
Lorentz, invariant. The difference between relativistic and non-relativistic strings is 
in the transverse perturbations. In an effective description, non-relativistic strings 
can be thought of as having a short wavelength cut-off on the string excitations. As 
a result, arbitrarily small- wavelength relativistic perturbations are not excited and 
this translates into a reduced curvature parameter normalisation /cq. The string can 
be thought of as a massive rigid rod, but with tension T equal to its mass per unit 
length n ^. In analogy to the relativistic case, where the overall normalisation was 



determined by comparison to a known analytic solution ||32|, k^ can be obtained in 
the non-relativistic case by comparison to a given model of non-relativistic string. In 
the general discussion below we will simply treat it as a free parameter and examine 
its effect on the network evolution. 



Equations (|6.24|) , ( |6.36| ) and (|7. 1| ) have been solved numerically for a range of 



parameters fco and W. This was done by rewriting equation ( |6.24|) in terms of the 
correlation length L = \/JjJp and then introducing a function 7(t) = L/t. Under the 
assumption L ~ 7?., the resulting equation for 7(t) together with ( |6.36| ) form a non- 
autonomous system of coupled ODE's, which can be integrated numerically. During 
matter or radiation domination, this system has an attractor solution in which both 



^Relativistic invariance in the longitudinal directions implies that the waves along the string 
travel at the speed of light c . Note the difference to the other notion of non-relativistic string 
with T < fi and longitudinal speed v < c. 
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Figure 1: Relativistic versus non-relativistic network evolution. Non-relativistic networks 
evolve to slower and much denser scaling configurations than their relativistic counterparts. 
Here, we have plotted the evolution of the normalised string density and rms velocity for 
a relativistic network and a non-relativistic one with V ~ 0.1. 



7(t) and v{t) tend to constant values (scaling). Here, we present numerical results 
for a radiation dominated universe. 

In Fig. [I| we plot the evolution of the string energy density and rms velocity 
for both a non-relativistic and a relativistic string network, that is, the solution of 
equations (|6.24| ), ( |6.36|) , (|7. 1|) in the former case and ( |6.5D , (|6.14|) , (|6.15|) in the latter. 
To highlight the effect of non-relativistic velocities, we have chosen a value of the 
parameter /cq which gives a scaling value of ~ 0.1 and taken W < V . We have 
also assumed that both networks have the same loop production efficiency parameter 
c and chosen the value c = 0.23, suggested by relativistic network simulations [p^] . 
As expected, the non-relativistic network has a much higher scaling string density 
compared to the non-relativistic one. Of course, non-interacting strings (c = 0) do 
not converge to a scaling solution. 

We now explore the dependence of non-relativistic string evolution on the pa- 
rameters ko and W. In Fig. |^ we plot the normalised string density pt^//i = 
and the rms string velocity V as functions of cosmic time t for different choices of 
producing string velocities < \^ < 1. We have assumed a constant value oiW <V , 
but below we will consider the effect of varying W also, allowing for W > V . It is 
apparent from Fig. ^ that the rms string velocities are controlled by the parameter 
/cq. Reducing /cq leads to smaller V , which in turn implies a higher string density, due 
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Figure 2: Evolution of normalised string density and rms velocity for a non-relativistic 
network with constant W < V for different choices of the parameter /cq. Reducing the 
value of ko results to lower rms string velocities, which in turns implies a slower rate of 
string interactions. This results in a dramatic enhancement of string network density. 



the reduced energy loss term. The fact that the scaling value of the rms velocity is 
not universal for non-relativistic strings, but instead depends on the parameter ko, is 
not surprising. In the relativistic case, there is a distinct upper speed limit c = 1 and 
the relativistic constraint implies that the rms velocities are smaller than, but not 
far off, 1/a/2 (see for example Ref. ^j)- On the other hand, in any non-relativistic 
theory velocities are unbounded. 

We then consider the impact of varying the parameter W. Looking at the first 
term of equation ( 6.24 ), which describes dilution due to cosmic expansion, one ob- 
serves that and appear with opposite signs, so a large W could counterbalance 
(or even reverse) the effect of V on this term. However, if both V,W <^ 1 they play 
no significant role in that term. Thus, one only needs to check the case W > V 
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Figure 3: Dependence of normalised string density on the parameter W for a network 
with V ~ 0.1. The plots correspond to W = 0,0.1 and 0.5 respectively. Increasing 
W does not significantly alter the scaling density until W becomes greater than V . For 
W = 0.5 = , the scaling is reduced by 10%, so it remains two orders of magnitude 
greater than that of relativistic strings. 



when y, W are not negligible. Fig. ^ shows the time evolution of p for a choice of fco 
leading to ~ 0.1, for the cases W = 0.1, 0.5. The first two figures show identical 
evolutions, even though in the second one ~ l^. In the third case, however, where 
= 25V'^ the effect of W counterbalances that of V in the dilution term of ( |6.24| ), 
resulting in an appreciable reduction of the string scaling density, at the 10% level. 
Since the most important impact of string velocities is through the loop production 
term of ( |6.24| ), the basic prediction of the model, which is a dramatic enhancement 
of the string scaling density (Fig. ffl), remains robust. 
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8. Discussion 



So far we have studied the dynamics and macroscopic evolution of non-relativistic 
strings in some generahty, without discussing any specific setup in which they could 
be relevant. However, non-relativistic string-like objects arise in several contexts and 
have been considered before in the literature. 



For example, Ref. |^ studied non-relativistic vortex-strings with motivations 

In that reference, 



from both cosmology [17| and condensed matter physics [42, 43 



the non-relativistic limit was taken at the level of the equations of motion by requiring 
small string velocities X"^ -C 1. Here, we have taken the non-relativistic limit at the 
level of the string action but this involved a rescaling procedure which corresponds 
to having both X"^ -C 1 and {dX/dC,Y ^ 1? where C, is the physical length along the 



string. The non-relativistic evolution model we have developed in section can be 
applied to the condensed matter context considered in |¥T] by introducing a friction 
term relevant to that situation. Adding this term and setting d/a = equation 
( |6.24|) , expressed in terms of the correlation length, reads: 

dL L o , , 



as in Ref. where id is the relevant damping length scale. The velocity evolu- 
tion equation is also modified by the addition of a friction term —id/L, again as in 
Ref. [^]. The system has a solution with L oc t^^^, which is actually observed exper- 
imentally for defects in condensed matter systems and liquid crystals ^ . 



In cosmology, slowly-evolving string networks have been invoked in order to 



obtain a negative equation of state Bucher & Spergel |^ have proposed a 

Solid Dark Matter (SDM) model, which could be realised in terms of a frustrated 
string or domain wall network. Rigidity and stability in this scenario have 
been studied in Ref. 



21 



More recently, a string network of the SDM kind was 
revived |22] in an attempt to explain the flat rotational curves and the TuUy-Fisher 
relation observed in galaxies, which were the main motivation for the development 
of MOND^" theories. The fundamental difficulty ^ with the SDM scenario is to 
explain how an essentially non-relativistic network can naturally arise from an initial 
tangle of (relativistic) Nambu-Goto strings like the ones beheved to be produced in 
cosmological phase transitions. Indeed, Hubble damping is inefficient at subhorizon 
scales [0 and there is no known mechanism efficient enough to damp the relativistic 
short-scale excitations on strings. These affect the equation of state through the 
velocity dependent term in equation (|5.3|) , leading tow > —1/3^^. Further, numerical 



^°For a recent review on the MOND scenario see Ref. p9| . 

"'^'^One could argue that the velocity which enters the equation of state is the coherent string 
velocity at the scale of the string correlation length rather than the rms short-scale velocity. While 
it is true that the coherent velocities are typically smaller, numerical simulations |I4| suggest Wcoh — 
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evidence is now accumulating supporting that scaling behaviour in field theory strings 
and domain walls is rather generic 0], so that frustrated networks seem hard to 
obtain. On the other hand the analysis we did in section |^ points towards a SDM 
picture for non-relativistic strings, where the above problems are not present. Here, 
string velocities can be arbitrarily small and, as we saw in section |^, network densities 
are dramatically enhanced so that strings could even dominate the universe before 
scaling is reached. 

Note that the procedure for obtaining the non-relativistic string action 
required at least one of the spatial directions to be compact. If the action (|3.11| ) is 
to be treated as a classical effective action this global property can be ignored, but 
if it is taken to describe a fundamental object, then the winding around a compact 
dimension is required at quantum level. The fact that a consistent non-relativistic 
string theory based on the action ( |3.11| ) can be constructed allows one to take 
the view that there is a fundamental winding string obeying this action. Then, a 
cosmological setup like that of sections and |^ can still be considered as long as 
the compactification radius is larger than the horizon. This possibility of having 
a universe with non-trivial topology is not observationally excluded. Cosmological 
observations constrain the local geometry as described by the metric to be nearly 
fiat [T^, but the global topology of spatial hypersurfaces need not be that of the 
covering space. Indeed, topological identifications under freely-acting subgroups of 
the isometry group are allowed, and the WMAP sky maps appear to be compatible 
with finite fiat topologies with fundamental domain significantly greater than the 
distance to the decoupling surface (see also [0]). 

One can therefore imagine a situation where fundamental non-relativistic strings 
are wound around 1-cycles in a non-simply-connected universe, in a setup analogous 



to that of the Brandenberger-Vafa scenario [^. If the compactification radius is 



larger than the horizon, as required by cosmological observations, a network of such 
wound strings behaves like an open string network. An analogous situation occurs 
in ordinary cosmic string simulations, where the network evolves in a periodic box 
and there is a class of long strings (determined mainly by initial conditions) which 
wind around the box. As the universe expands these strings tend to straighten out 
and behave essentially non-relativistically These strings are usually discarded 
as artifacts of the periodicity of the box, but in a universe of compact topology, such 
configurations can play a physical role. 

Finally, in theories with compact extra dimensions one has the possibility of non- 
relativistic strings winding 1-cycles in the internal space. Analogous (but relativistic) 



objects have been considered in the context of brane infiation [55, Bq, pTj, which are 



topologically trapped and behave like monopoles. Although the copious production 



0.15 so one still expects significant departure from w — —1/3. Furthermore, small-scale structure 



has the effect of 'renormalising' the string mass per unit length |17, uM and string tension so that 



equation (5.4) should be used instead of (5.3). This also increases the value of w. 
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of such objects in the early universe is inconsistent with the existence of an early 
radiation era, there are regions in parameter space where they are allowed and in 
some cases can provide candidates for dark matter. The situation of non-relativistic 
strings wrapping an internal dimension is qualitatively similar, but the corresponding 
energy spectrum is different than in the relativistic case. 

The outstanding question arising from the above is to what extent such non- 
relativistic strings are 'natural' or 'generic' objects in cosmology. Even though non- 
relativistic strings exist in some part of the moduli space of string theory, there is at 
present no mechanism which produces them in a cosmological setup. Nevertheless, it 
is clear that the non-relativistic string action and the VOS model developed here are 
applicable at least as effective descriptions of cosmic- and vortex-strings in certain 
situations. Indeed, the action we have considered is the only sensible non-relativistic 
limit, having T = /i, of the standard Nambu-Goto action, and is precisely the ac- 
tion one obtains when considering the low energy dynamics of topological defects in 
field theory. The macroscopic NRVOS model based on this action, provides a semi- 
analytic tool for the study of the cosmological evolution of non-relativistic strings. 
Possible situations of cosmological interest involving non-relativistic strings include 
strings in de Sitter space, Sohd Dark Matter, wound strings, etc, as discussed above. 
Further, in a condensed matter application we have noted that our model reproduces 
the correct scaling law, as experimentally observed. 

It would be interesting to go one step further and perform numerical simulations 
of string network evolution based on the non-relativistic string action presented here. 
The comparison of macroscopic string evolution and small-scale structure to the 
relativistic case could provide an independent means of probing the effect of small- 
scale structure on string networks, which is an area of current interest and active 
research. 
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A. Non-Relativist ic String from Semiclassical Approxima- 
tion 



In this section we derive the action (|3.11|) as a semiclassical expansion around the 



vacuum solution. Non-Relativistic D-brane actions on AdSs x have been recently 
constructed with this method, by considering the Dirac-Born-Infeld (DBI) action 



and expanding around a classical solution. Here, we apply the method to the case 
of the Nambu-Goto action in an expanding FLRW background. 
We start from the Nambu-Goto action 



Sng = -TJ V^d'a (A.l) 

and write the induced metric •jij = GuNdiX^^ djX^ as 

lij = efefvMN , (A.2) 

where ef^ = a(x°)(9jX^^ is the (worldsheet induced) vielbein for the FLRW metric 
). We now consider the vacuum field configuration: 

r , M = 

x^^ = { Xa ,M = 1 (A.3) 
, M = ae (2, 



which is a solution of the equations of motion ( ^2.21 ). The only non-trivial vielbeins 
evaluated on the solution are the longitudinal ones 

e^i: = a(x°)9,X^ /i = 0,l, (A.4) 

the transverse ones Cg'^j being zero. The induced metric on the static solution is then 

7oy = eoie^jT]MN = eo^ej^-r]^^ . (A.5) 

Next we introduce transverse fluctuations around this solution, namely 

= X'^ia') . (A.6) 

The induced metric can then be written as the sum of the static solution plus a piece 
quadratic in the fluctuations: 

lij = loij + l2ij , (A. 7) 

where 

^0^, = a{xydia.g{-{^rxy, (9^x^)2) =a(a;°)Miag(-l,A2) (A.8) 
l2^J = a{xy^,X%XX,. (A.9) 
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Expanding ^ we have: 



\/-det7 = A/-det(7o + 72) = V -det[7o(l + 7o S2) 

1 



- V -det7o \ \ + -7o^72ij + • 

= A/-det7o + ^A/-det7o7o^72y + . . . . (A.IO) 

For the zero-order (in the fluctuations) piece we can write: 

V-det7o (iV = det(eo)f(iV = ^£^,,60%'' , (A. 11) 

where our conventions are such that cqi = 1. This can be therefore cancelled by 
choosing a closed B^^ field as in section ^ We are left with the second-order piece, 
which corresponds to the non-relativistic string action ( p. Ill ) we obtained in the last 



section. Indeed, from equations ( |A.1| ), ( [A .91 ) and ( |A.1(J| ) we get 



S2 = -^j a{x'Y^/^d^o7o'SabdiX''djX''d^a = j v/=d^7o'Ga6(x°)9,X'^6',X^ciV . 

(A.12) 

B. Hamiltonian Formulation 

Here we discuss the Hamiltonian formulation of the Nambu-Goto and non-relativistic 
strings. 

B.l Relativistic String 

We first consider the relativistic Nambu-Goto string. The Hessian of the Lagrangian 
( |2.1| ) has two null eigenvalues and as a result the canonical variables x^, Pm = ^q±m 
satisfy two primary constraints, namely 

PMpLG'^^(x') + Tx'^'x^'GMNix') = (B.l) 
PMx'^' = 0, (B.2) 



which are first-class. From these the Dirac Hamiltonian can be constructed 



H = / Hda 



da 



(B.3) 

where the Lagrange multipliers f,g are arbitrary functions on the worldsheet. The 
Poisson Brackets for x^ and pm are 

{X'\a), Pn{cj')} = S{a - a')5^ (B.4) 
{X^^(a),X^(a')} = (B.5) 
{PM{cT),PN{(r')} = (B.6) 
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and, choosing h = 0, the equations of motion read: 



X 



M 



{x^', H} = fT-^pN{a)G 



PM = W, H} = T + {T-'pnPaG'^H^') - Tx^'x'^GM^')) S% 

(B.7) 

Thus, for the spacehke fields x^ , I = 1, . . . ,D, we have: 

= / {fx'^)'6i - 2^fT''pja{x')-'6'' + fT-'pja{x'r'6^ 
^ a 

In the temporal gauge i'' = 1 we have from equation (|B.7|) that 

/ = -Ta(x°) V' 

and so 

Pi = Ta{x^Y f~^x^ = —pox^ . 
Then equation (|B.8D becomes 

x' = T'a{x')%' {p^'x")' - poPo'x' . 

Now po can be found from the constraint (|B.1[ ): 



(B.8) 

(B.9) 
(B.IO) 

fB.ll) 



Po = Ta{x 



0\2 



X 



1 — X 



12 



Ta{xye 



(B.12) 



In view of our previous results for e in the Lagrangian formulation (see equation 
(P^.6|)) we see that we are going to recover the equation of motion for x^ . Indeed, 
from ( [B.Tp we have: 



Po % = -aaTpQ ^ T ^a(a;°) ^ {plx^'^ - pi) - Ta{ 



2- (1 - x'') 



where we have used ( |B.12 ). Thus, ( B.ll ) becomes: 

{e-'x")' -2^{l~x'^) 



x' = ^^-'-"^ 



(B.13) 



(B.14) 



as in (|2.6|) . Finally, the constraint ( p.2|) in this gauge becomes x^x'^ = 0, as before. 
B.2 Non-Relativistic String 

We now turn to the Hamiltonian formulation of the non-relativistic string. The 
Hamiltonian density can be constructed from the constraints ( |4.1[ - |4.2| ) introducing 
arbitrary functions /(a, r) and h{a,T) as Lagrange multipliers: 



n = f\p,e'^%y'' + - 



1 rPnB 



Tr 



C'^ix'') + 







(B.15) 
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The canonical variables satisfy the following Poisson Brackets: 



N 



{X^^(a),X^(aO} = 
{PM(a),PAr(a')} = 0. 

Then, choosing h = the transverse equations of motion are: 



(B.16) 
(B.17) 
(B.18) 



X'^ = {X^ H} = fT^^P,{a)G''\x'') 
Pa = {Pa, H] = To {fX'\a))'Gat{x') 



(B.19) 



while the longitudinal ones read: 



(B.20) 

Now choose the gauge ( [4.5| ) by setting x^ = 1 and x'^ = A. Equations ( |B.2(J| ) then 
imply 

/ = -A-^ (B.21) 



and the transverse equations of motion ( [B.19|) give 



X'^ = -(ATo)-^ (^PhG''\x'^) + PbG^^x^' 



recovering equation (|4.7| ). 

From equations (p.20| ) we have: 



Po = A-i 



/ 0- ( PaPb^ab(^0 
Pi 



a \ Tf) 



G'^'ix'') - ToX"'X"'Gab{x°] 



a 



^G''\x^] 

Tn 



(B.22) 



(B.23) 
(B.24) 



using the constraint ( [4.2|) . Since po is given by equation ( [4 .11) and Pa = — AToX^ 
(from (p.20| )) we recover equation ([4.13|) . 



C. Energy-Momentum tensor of Non-Relativistic Action 

In order to obtain the energy-momentum tensor of the non-relativistic string, we 
vary the action (|3.11|) with respect to the background metric Gmn- 



T 



MN 



ss 



V— det^MAT 6G 



ss 

SGu,„ 



MN 



V-detG 



MN 





ss 

SG„h 



(C.l) 
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The transverse part reads 



To 



V-detGMAf 
whereas the longitudinal part is: 

"1 



J £a^/^d^g''diX''djX''6^''+^\y^' - x^\a')) (C.2) 



^J -detG MN J 



2 



-g'"'d,X''dmX''g^^d,X'dnX''Gab] 5'~''+'\y^' - x^' {a')) (C.3) 
In particular, the 00 component in the gauge ( |4.5| ) becomes 

(C.4) 

where, having integrated out S{ri — t), K runs from 1 to D. Note that the ex- 
plicit scale-factor dependence of the integrand cancels, because of the presence of 
y^^g~^g~^G in equation (|C.3| ). The first factor scales like airjY, the next two fac- 
tors as a{ri)~'^ each, and the last factor as a(?7)^, giving a scale-factor-independent 
result. Of course T°° still depends on time through the time dependence of the fields 

x^ 

The string energy can be defined as in section |], by considering a spacial hyper- 
surface ?7 = const, with normal (co)vector = (a(?7), 0), and integrating the energy 
density hmtinT'^^^ = a{riyT^^ = Tq over the relevant D- volume: 

E{r]) = -J VhnMnNT^''d''Y 

= a{7]) ^Jda (X^'X' + A^'X'^X'^) 6ab • (C.5) 
Similarly, the ab components of the energy-momentum tensor in the gauge ( |4.5| ) 



are 



T^\v,y'') = ^^^, j\To{-X'^X' + \-'X'^X")6^''\y'' -x^{a,r^))da. (C.6) 
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